Analytical and numerical solution of coupled KdV-MKdV system 
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Abstract 

The matrix 2x2 spectral differential equation of the second order is considered on x in (— oo, +00). We 
establish elementary Darboux transformations covariance of the problem and analyze its combinations. 
We select a second covariant equation to form Lax pair of a coupled KdV-MKdV system. The 
sequence of the elementary Darboux transformations of the zero-potential seed produce two-parameter 
solution for the coupled KdV-MKdV system with reductions. We show effects of parameters on the 
resulting solutions (reality, singularity). A numerical method for general coupled KdV-MKdV system 
is introduced. The method is based on a difference scheme for Cauchy problems for arbitrary number 
of equations with constants coefficients. We analyze stability and prove the convergence of the scheme 
which is also tested by numerical simulation of the explicit solutions. 
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•^1 Introduction 

• 1— 1 . 

Xxhere are two complementary approaches to integrable systems: analytical and numerical ones to be devel- 
odoped. Even most profound analytical 1ST method cannot give explicit solution of general Cauchy problem 
while numeric can, but is rather compulsory in use: calculations could need powerful computers. May be 
most transparent of analytical methods are based on algebraic structures associated with a problem. To 
such structures belongs Darboux transformations covariance of Lax representation of nonlinear equations 
that yields a powerful tool for explicit solutions production. We investigate applications of special kind of 
such discrete symmetry - to be called elementary ones ||10|| . Its elementarity simply means that a product 
of such transformations generate the standard one (TO, Here we study combinations of such transforms 



that do not coincide with binary ones []T2| and hence are not so known. 

The main ideas of numerical integration of such integrable systems go up to the famous properties of 
the equations as the Lax pair and infinite series of conservation laws existence [[!]]. From a point of view 
of general theory of such systems some hopes are concerned with a development of the finite-difference or 
other approximations of the systems. Namely if one could prove a convergence and stability theorems for 
such difference systems (existence of solutions is implied), a way to existence and uniqueness of solutions 
is opened ||. 

The coupled KdV-MKdV system arises in many problems of mathematical physics. Some integrable 
systems are associated with a polynomial spectral problem and have Virasoro symmetry algebras are 
considered |[L7| . A dispersive system describing a vector multiplet interacting with the KdV field is a 
member of a bi-Hamiltonian integrable hierarchy ||. Recently a multisymplectic numerical twelve points 
scheme was produced. This scheme is equivalent to the multisymplectic Preissmann scheme and is applied 



to solitary waves over long time interval [O]. The coupled KdV-MKdV system is also connected to other 
physical applications |TI ]. 

The general system we consider in this work have the following form 



+ E {& m e k x + g n m 2 k (e m m + + O m C + tf fc <H* fc )a + 4C = o, (i) 



rn,k 



where n, m, k =1,2, ...N are the dependent variables numbers. Nonlinear coefficients are g^ k , I =1,2, ..5 
and d n are dispersion coefficients. 

In particular, for the system under consideration (N = 3) the variables 1 , 2 , 9 3 are denoted by f, u, v 



to have the integrable system [10 



ft + -fxxx + -^{ u f)x — -fxf 2 — 

u t - -u xxx - ^u x u + 3vv x + ^u x f 2 - ^(f x v) x = 

1 3 3 3 3 

v t + -Vxxx + -jV x ii - - {vf) x + -u xx f + -u x f x = (2) 



The Lax pair is given in flnj. The system exhibits two integrable reductions having explicit solutions, 



Hirota-Satsuma @, [l|] and a two components KdV-MKdV system |10|, [121 . Krishnan || showed that a 



generalized KdV-MKdV system have solitary wave solutions and investigate the effects of increasing the 
nonlinearity of one variable on the existence of solitary waves. Some form of KdV-MKdV system have 
explicit solutions in terms of Jacobi elliptic functions |2j . 

In this work we present explicit solutions for a system of three equations (0) that have not been specified 
TP], O. We study this two-parameter explicit solutions and show effects of choosing these parameters 



m 



on the solutions. We demonstrate the use of two arbitrary elementary DTs [10| and its special choice that 
holds a hereditary of the reduction to built explicit solutions to the KdV-MKdV system (^). 

Also we modify a numerical method |7], |3[] for solution of system ([!]). It is a difference scheme for 
Cauchy problems for arbitrary number of equations with constants coefficients. The scheme preserves two 
conservation laws for the KdV type equations and the order of error of the difference formulas is improved 
0, 15]. The convergence is proved and stability is analyzed giving the conditions taken in account in 



choosing time and space step sizes |[6|, . 

The present work is organized as follows. Section 2 introduces the matrix spectral equation of the 
second order with 2x2 matrix coefficients and two elementary DTs. We select the second equation of 
the Lax pair and derive the compatibility conditions. The product of these two transformations yields 



the standard DT |10|, [13| . Section 3 illustrates how the, first, elementary DT is used to produce solution 
to the KdV equation as well as the general evolution equation generated by the compatibility conditions 
of the Lax pair. Explicit solutions are introduced for the case of zero initial potentials of the matrix 
problem. In section 4 we consider a reduction constraints on the potential of the matrix spectral equation. 
This reduction gives an automorphism that relates two pairs of solution of the spectral equation for two 
spectral parameters. We use this results in the compound elementary DTs to produce an explicit solution 
to a coupled KdV-MKdV system that results from the compatibility conditions of Lax pair under this 
reduction. The effects of these parameters on the solution (reality, singularity) is analyzed. Section 5 
introduces a numerical method for solving coupled KdV-MKdV system (U). We produce a difference 
scheme for a Cauchy problem with initial condition rapidly decreasing at both infinities. The main steps 
of the scheme convergence and stability analysis is shown while the details are explained in appendix A 
and B. The scheme is tested by applying it to integrable coupled KdV-MKdV system and the numerical 
results are compared with explicit formulas obtained in section 4. 



2 Lax pair spectral equations and the elementary DTs 

Consider a matrix spectral equation of the second order with spectral parameter A and 2x2 matrix 
coefficients. 

^ xx + FV X + U* = \a 3 ^ (3) 

where the vector \1> = \l/ 2 ) T and the matrix potentials are U = {uij}, F = {fy, fa — 0}, i — 1, 2 while 
o"3 = diag(l, —1) is the Pauli matrix. 

For equation @ we perform two elementary Darboux transforms [ 10| . The first one is 
#i = + en^i + £12^2 , en = - ( ip lx + \j\2f2 ) /fi , e i2 = /12/2, 



2* 

^2 = ^2 + £21 ^1 , £21 = - <£>2/Vl> 

£1 = (<9x + en) y? 3 + e 12 <£>4 , ^2 = ^4 + 621^3- (4) 

where (</?i, v?2) T and (</? 3 , v 9 4) T are two solutions of (|3]) corresponding to different spectral parameters. 

Substituting the above expressions for ^ 2 into (0) and collecting the coefficients of ^1, ^! 2 and their 
derivatives we obtain the expressions for new potentials as 

/12 = W12 + /12 en, 
/21 = — 2 e 2 i, 
itxi = "11 — 2eiia; — f l2 €21 — /21 ei2, 
Wl2 = Ul2x - £l2xx + en«i2 - e i2 (u u + U 22 

u 2\ = /21 — 2e21x — f 2 \^U, 

U22 = u 2 2 - e 2 iwi2 - «2iei2 - / 2 i e i2x- (5) 

The second elementary DT is performed after the first one and can be obtained by reversing the indices 
1— > 2 and 2— > 1 to get, for example, the following potentials 

/21 = «21 + /2ie22, e 22 = - ($2x + ^2l£lJ /<f>2, 

U22 = u 2 2 - 2e 22a; - / 2 iei 2 - /i 2 e 2 i, e 12 = - £i/£ 2 > e 2 i = f 21 / 2. (6) 
The spectral equation (0) is considered as the first equation of the Lax pair, take the second as 

% = * xxx + BV a + CV, (7) 

where B = \diagU + \F X + \F 2 

and C = \U X - \diagU x - |(/i 2 m 2 i + /21W12K + K/i^^i - fuhi^a + f («ii - u 22 )a 3 F. 

Equation (0) is also covariant under transformations (H). The compatibility conditions have the 
following form 

F t - F 3x . + S 2 , - W 2x + 2C X + FB X - a 3 Ba 3 F x + [75 

-a 3 Ba 3 U + FC- a 3 Ca 3 F = 0, 
U t - U 3x + C 2iE + UC- a 3 Ca 3 U + FC X - o 3 Bo 3 U x = 0. (8) 

and the transformations (PS), (H) determine a discrete symmetry of (0) 



3 Solution of two coupled KdV-MKdV equations and KdV 
equation via the first elementary DT 

For a spectral parameter A and a seed potential F, U we obtain the solutions (fi,(f2 to the pair (|3|), (0). 

Then performing the first elementary DT to obtain the new potentials F, U which are solutions to the 
system (§). For the case of zero seed potential the solutions, <pi and <p 2 of the system (f|), (|7|) have the 
form 



ax+a^t i — (ax+a 3 t) 

tp i — c\e + c 2 e y ', 



¥>2 



die iax+(*a)H + ^-(iax+MH) _ ^ 



where c±, c 2 , d\, d 2 are arbitrary constants, a = \/X and i is the imaginary unit. System (|2|) reduced (for 
the only nonzero elements ) to the following 

J2Xt ~r -J21xxx ~r ~J21 v 'llx — —-U11U21, 

1 3 -n 

Ullt — —Un xxx — — U\\U\\ X — U, 
«21t + ~V>21xxx + ^U 2 lU Ux + ^U 2 lxUn = - U n f 2 l xx + jWn/21- (10) 

where fi 2 = 0, U12 = 0, u 22 = and tildes are omitted for simplicity. This system with explicit solution 
obtained from d|), @ as 

f 21 = (2e^- l) < aH+x ) [d 2 e 2iaH + d x e 2i ™\\ I (c 2 + Cl e 2 < a2t+X < 
u u = (8a 2 c lC2 e 2 < a2t+x )) / (c 2 + Cl e 2 < a2t+X ) 
U21 = (-2iae {1 - i] < a2t+x ) [d 2 e 2iaH - d x e 2iax ^ / (c 2 + Cl e 2 < a2t+X ^ . (11) 

where ci, c 2 , c?i, d 2 are arbitrary constants and a = y/\. 

The second equation in ( |I0| ) is the KdV equation while the remaining are a two components coupled 
KdV-MKdV system that was solved by elementary DT. 

4 Solution of three coupled KdV-MKdV equations via the com- 
pound elementary DTs 



2 



Existence of different kinds of automorphism causes special constraints [10]. Multiplying (Q) by ci = 
(l 0) ^° ^ ave 

(ri^xx + (TiFV x + a 1 U^> = \a 1( T 3 ty (12) 
but o"iO"3 = —0301 and consider the conditions a±F = Fa\ and o\U — Ua\ that means 

f\2 = hi = /, un = u 22 = u, U12 = U21 = v. (13) 

So ( |12"D becomes 

( a ^)xx + F (^1*)* + U ( CT i^) = (^1^) 



The above automorphism \1/ (A) <— cr^ (—A) relates two pair of solutions (fi,<f2) and (<^ 3 ,<^ 4 ) of (|3|) 
corresponding to different values of spectral parameter A, —A as 



<p±(-\)J ai \ip 2 (\)J \fi(X) 

Using this result in the elementary DTs (|j), (|5|) and (|6|) to obtain the expressions for the new potentials 
f, u, v. In the case of zero initial potentials these new potentials have the following forms 



M 2 - (v> 



2) 



_ , (<Pi)x - (<pZ)x \ +2 f <piMx -^W 



M 2 -M 2 J' \ M 2 -(<p 



J 2 



2 



fo>i) 2 -M 2 Jx (M 2 -M 2 Y 

where (pi, (f 2 are as in (|9|) with ci, c 2 , d 2 are arbitrary constants and a = y/X. 

The above expressions are solutions of system (§) that reduced under the reduction conditions (fl3|) to 
system (Q). 

The choice of the arbitrary constants (ci, C2, di, d 2 ) affects on the behavior of the solution in formula 
(|j). For example choosing equal constants ci = c 2 = c?i = d 2 = 0.5 (we choose the value to be 0.5 to 
simplify the resulting formula but the idea valid for any value) the solutions have the form 

/ = 2a(sinrjicoshrj 2 — cosrjisinhrj 2 ) / \cosh 2 r] 2 — cos 2 rji), 
u = 2a 2 (sinr]iCOshr] 2 + cosr\\sinhr) 2 ) 2 j \cosh 2 r] 2 — cos 2 r)i) 2 , 
v = 2a 2 (cos3r]icoshr] 2 — 2sinr]isinhr] 2 (cos2r]i + cosh2i] 2 + 2) — cosr)icosh3r] 2 ) 

j (cosh 2 i] 2 — cos 2 r]i) 2 . (15) 

where 771 = a 3 t — ax, r] 2 = a 3 t + ax, a = y/X is real. 



We see that the above expression ( |15|) is singular at r\ 2 — 0, r\\ — rvK,n = 0, 1,2, .... Hence we have 
singularity at (x = *g,t = 

To obtain continuous solutions we can choose c\ = c 2 ,d\ = d 2 = r.ci, r is real constant. We again 
choose c\ = 0.5 following the previous concept. So (|J) have the form 

/ = 2ar (cosh r\ 2 sin 771 — cos r\\ sinh 772) / (cosh 2 r\ 2 — r 2 cos 2 771) , 

u = a 2 (l — r 4 — r 4 cos 2i}\ + cosh 2772 + r 2 sin 2r\\ sinh 2772) / (cosh 2 r\ 2 — r 2 cos 2 771) 2 , 

v = 2a 2 r(((— 7 + 6r 2 + 2r 2 cos 2?7i) cos 77! cosh 772 — cost/! cosh 3772) 
—2(1 + r 2 + r 2 cos 2r/i + cosh 2i] 2 ) sin 771 sinh r] 2 ))/(— 1 + r 2 + r 2 cos2r/i — cosh 2??2) 2 (16) 

where a, 771, 772 as in fllED 



Choosing this parameter (r) to be r < 1 gives real nonsingular solutions. The above formula (f|) is 
built from elliptic and periodic functions so it does not preserve its symmetry but its localized as shown 
in figures (l.a,b) below. 



Fig.(l.a) Non-singular solutions, f, u and v (r=0.5) , a=2, t=0. 




Fig.(l.b) Propagation of solutions, f, u and v (r=0.5) , a=2, t=l. 
Fig.(l) The solutions in does not preserve its symmetry but its localized. 

Choosing the parameter r to be r > 1 in formula (H) gives singular solutions as shown figure (2) below. 




Fig. (2) Singular solutions, f, u and v (r=2), a=2, t=0. 

Moreover the choice of these arbitrary constants (c±, C2, d±, d 2 ) as well as the spectral parameter A affects 
the reality of the resulting solution. As example for A = —2im 2 , m is real and choosing c\ = c 2 = d\ = 
d 2 = 0.5, we get real solution 

/ = m(cos2(isinh( 2 — sinh( 2 — sindcosh2( 2 + sind) / (0.25cosh2( 2 — 0.25) (1 — cos2d) 

where (i = 2mx + 4m 3 t, ( 2 = 2mx — Am 3 t, while choosing C\ — c 2 = 1, di = d 2 = 2 give the following 
complex solution 

/ = m * (—8(—5sinh( 2 cos2(i + 5sinh( 2 + 5sin(icosh2( 2 — 5sin(i) — Si(6sinh( 2 
+3sin2(icosh( 2 + 3cos(isinh2( 2 + 6sin(i)) / (17 cosh2( 2 + 10 + 36cosdcosh( 2 

-8cos2( 1 cosh2( 2 + 17cos2Ci) 



5 The numerical method 



5.1 The difference scheme 



For the coupled KdV-MKdV system ([I]) we introduce a numerical (finite difference) method of solution 
0. |l5| . This scheme is valid for arbitrary number of equations with constants coefficients and of the form 



an,j+l nn,j nk,j _ nk,j nk,j nk,j 

°i v i , V~~V n,l Q m,j u i+l . n,2 / Q m,j\Z u i+l 

+ 2J<9^i 2h + 9m ^ 1 > 2h 



m,k 

nm,j _ nm,j nk,j _ nk,j nk,j _ cynkj , nk,j 

, n,3 u i+l u i-l u i+l u i-l , n,4nm,j °i+l ~r 

+9 m , k 2h 2h +9 m ,k i 2h 

nk,j _ nk,j nn,j _ nnnj , nnn,j _ nn,j 

,_n,5 gm,j a k,j u i+l , j °i+2 zf7 i+l T zc [i-1 °i-2 _ n /-. 7 \ 

+9 m , k Vi f< 2h )+dn ^ -U [L() 

where i and j are the discrete space and time respectively. The time step is denoted by r while h denotes 
spatial step. 



5.2 Stability analysis of the scheme 



We prove stability with respect to small perturbations of initial conditions |16, [|. It is the boundness of 
the discrete solution with respect to small perturbation of the initial data. We give here the main steps 
while the details are presented in Appendix A. We can write 



<:W 



nj'+l 



T?; j+1 d9^ = T^T^de^- 1 = IT (7£) r d0?° 



where d6™' 3+1 is perturbations of the discrete solution, dd^'° small perturbation of the initial data and 
T" , ' J+ is a different iable operator. Stability required the boundedness of II (T^Y i.e ||T r || is bounded. 

' r ' 

We found that 



a(r, h) = 2 max 

Ln,m,k 



n,l 
9 m ,k 



max ( \9 rr . 

i,m,k 



xA 



k,j 
x,i 



+ r[ max 

l,n,m,k 



n,l 
9 m ,k 



max | 0%) 

i,m,k V / 



+ — max 

II l,n,m,k 



n,l 
9m,k 



max ( \6T' 3 

i,m,k 



+ — max|d n |] 2 

tl° n 



The scheme is stable if a(r, h) < const. We have here a conditional stability. That is we require that 
r — > more faster than h — > . Namely we need 



r < (constant) . ti 



5.3 Convergence proof for the scheme 

We prove that the solution of ( |17D converges to the solution of ([[]) if the exact solution is continuously 
differentiable one ]TB|, |S|. We introduce here the main points for the scheme convergence and give the 
details in Appendix B. 

Q\ is the difference solution of (|l7|) , u? is the exact solution. Hence the error v] is given by v\ = d\— u\ . 
Introducing L 2 norm defined by ||V^|| = £E (v? J ) h 

\ i n 

The scheme converges when the norm of that error ||V J '|| — > as (r, h — > 0) 
We found that IIV^ 1 !! < P(M) 0(r + h 2 ), where P(M) is a polynomial in the bounded constant 
M = t ea JJrZ\ an d a as in (|18|). Hence the convergence proved. 



i 



5.4 Numerical calculations and test 

The coupled KdV-MKdV system (Q) is solved numerically using scheme (p!7|) with initial condition from 
(|) at t — and the results are compared with the explicit formulas (|J). The percentage errors are shown 
in the following plots. 
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Fig. 3 percentage errors of the numerical solutions relative to the explicit solutions. 



The results of the test confirms the validity of the numerical scheme we propose. It also illustrates the 
errors of evaluation that could be estimated by the resulting inequalities of the scheme convergence proof. 



6 Conclusion 

Darboux transformations covariance of Lax representation of nonlinear equations is a powerful tool for 
explicit solutions production. Here we investigate applications of special kind of such discrete symmetry - 
to be called elementary ones. We use these elementary DT to produce explicit solutions for coupled KdV- 



MKdV system. The iteration of DT can be formulated in form of determinants representations |10|, 12] . 
A numerical method for general coupled KdV-MKdV system is introduced. It is a difference scheme for 
Cauchy problems for arbitrary number of equations with constants coefficients. We analyze stability and 
prove the convergence of the scheme. The scheme keeps two conservation laws chosen in analogy with KdV 
type equations. Analyzing stability and proving the convergence beside comparing the numerical results 
with explicit formulas allow us to use the numerical scheme to systems with arbitrary coefficients that is 
presumably non-integrable. Obviously the coupled KdV systems are successfully treated by our scheme ||. 

Acknowledgment We thank S. B. Kshevetskii for useful discussions about numerical scheme for 
the problem under consideration. 
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7 Appendix A 

Stability analysis of the scheme 

We prove stability with respect to small perturbations (because we consider nonlinear equations) of initial 
conditions. Strictly speaking it is the boundness of the discrete solution in terms of small perturbation of 
the initial data. Consider the differential 

T n ;j +i = { 9 0?j+i/d0nj} , d ffnJ = 0£J #^ 2 }Vd define the norm ||c^'|| = ( EE W ,j f h 

V r n 

We can write d9^ +1 = T^ j+l d9^ = T^ 1+1 T n ; j d9^~ l = IT (TfY d9?>° 

where is perturbations of the discrete solution, d9^° small perturbation of the initial data. 

Stability required the boundedness of Id (T™ r ) r i.e ||T r || is bounded. We calculate T from ( |T7| ) as follow 



-inj'+l 



m,k 



n,l 
■ 9m,k 

~2h 



n,2 
9 m ,k 



n,3 
9 m ,k 



2h 



2h 



n,4 
9m, k 



n,5 
9m,k 

' 2h 



2h 

m,j a k,j 



rd r . 



2 k J 



k,j I nk,j 



2h 3 



[8i+2,r ~ 2<5i+i. r + 25j_i )? . — 25j_2,r 



(19) 



Rewriting (BJ) in terms of identity (E), symmetric (S) and anti-symmetric (A) matrices 

o k x ' J 



\S^ + || < t max 

l,n,m,k 



9m\k 



ma ?(ICi J| 

i,m,k x 



\\A^\\ < I max 

Ln,m,k 



n,l 
9m,k 



max ( 1 9T' 3 

i,m,k 



x,i 
) k ,j 



+ ^max \d n \ 



where 91 



t + 1 i—l 

2h 



n,m,k = 1,2, ..N, I = 1,2, ...5. 



\( Tj+ iy Tj+ i\\ = \\( E _ A j+i + ( E + A j+i + si+i) 



<l + 2 ||5i+i|| + (||^ +1 || + ||5 , i+ 1 ||) 2 

< e a(r,fe)r 



a(r, h) = 2 max 



n,l 
9m,k 



max ( 1 9™f 1 

i,m,k 



r[ max 

Ln.m.k 



n,l 
9m,k 



max ( I 



i,m,k 



+ — max 

/i l,n,m,k 



n,l 
9m, k 



max ( | 



+ —max\d n \} 2 

tl n 



We have here a conditional stability. That is we require that r 
need 



more faster than h — ► 0. Namely we 



r < (constant) . /i 6 



8 Appendix B 

The scheme convergence 

We prove the convergence by proving that the norm of the error (between the difference solution and the 
exact solution) vanishes as the mesh is refined. Let 6{ the difference solution of (0), u^ the exact solution. 
The error vj is given by v\ = 9j — u\ . 

The scheme converges when the norm ||V 3 '|| — > as (r, h — > 0) where the norm is defined as \\V^\\ = 

(ee {<' j ) 2 hj 2 

substitute in ( |T7| ) by 9\ = v? + u\ keeping in mind that for Q\ equation (|TT| ) is O (r + h 2 ) and using 
the operator T defined by 



k.j k,i k.j k.j k,j k.j 

flj 7J 7/ 7/ 7 J 7? 

n,7 r\ V n,l / m,j u i+l i—1 , m,i : i+1 i— 1 \ , n,2 // m,i\2 i+1 i—1 

Vi ~ r[}^(g m \k( u i Vi 2h ) W( u i ) 2ft 

fc,j m,j m,j k,j m,j m,j k,j k,j 

, / m,j\2 U i+l ~ U i-1 \ , n,3 ^i+l ~ %-l V i+1 ~ ^i-1 . ~ V i-1 U i+1 ~ %-l ^ 



2/i ' * m>fcV 2h 2h 2h 2h 

, n,4 / m,j W i+1 — ^ V i > V i-1 i m,j U i+l~^ U i ' U i-1\ . n,5 / m,7 

+S'm.fcK SI + ^ SI ) + 2 9 m \k( u i v , 



k,j k,j k,j k,j n,j c\ n,j . r> n,j n,j 

V i+1 ~ V i-1 , m,j k,j U i+l ~ U i-1 \\ , j ^+2 ~ + ~ V i-2 ] 

2h +Vl Ui 2h ))+ n 2h? J 



j 



r 



So we obtain 

<•:"•'•' <•;:••' • rd, (20) 

r 

k.j k,j k.j k.j m.j m.j k.j k.j 

i rn,j n >l m )0 v i-l I / m,'j \2 V i+l~ V i-l I n >3 v i+l~ v i-l v i+l~ v i-l 

Wnere Jm,k,i — 2^9m,k V i 2h ' 9 m ,k\ V i ) 2h ' 9m,k 2h 2h 

m,k 

k.i r, k.i i k.j , fc , 7 fc , ? 

n,4 m,7 ^+1-2^ . n,5 m.j ftj^i+i - ' v i-l , r\ / i !,2\ 



\\p 



v(./:i,f /> 



i 71,1 1 

<- l£mi* I max 

- P72 

Ln,5 I 



1/2 



I "-,2 



I n i3 



n,4 
Qm,k | 



< 



h 2 

I I 

9ml 



' + ™^ ll^'f + li^J-* ll^'f + ™^ ll^ll 

+ O (r + /i 2 ) 
+ O (r + /i 2 ) , 



n,l 
9 m .k 



n,rrt,fe,l ' 



Using Schwartz inequality so ( |20| ) becomes 

ll^ i+1 ||<ll^ +1 || + t\\P\\ 

< || T i|| ||^i-i|| +r(T i+i|| / ,--i|| + 11^11) 

< e aTj ||y°|| + r (e^^- 1 ) ||/°|| + e aT ( J '" 2 ) H/ 1 !) + .. 



< e aTj \\V°\\ +M 



\\f j \\) 



9 m ,k 



\\V j+1 \\ +M 0(r + h 2 ),M = r 



Using ||U°|| = 0, the above inequality has the solution 

< P (M) O (r + h 2 ) , P(M) is a polynomial in M. 
Since M is bounded then ||V- ? ' +1 || — > as r, ft, — > and the convergence proved. 



